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Abstract. In this article it is shown that in a classical equilibrium canonical ensemble 
of molecules with s-body interaction full Gibbs distribution can be uniquely expressed 
in terms of a reduced s-particle distribution function. This means that whenever a 
number of particles N and a volume V are fixed the reduced s-particle distribution 
function contains as much information about the equilibrium system as the whole 
canonical Gibbs distribution. The latter is represented as an absolutely convergent 
power series relative to the reduced s-particle distribution function. As an example a 
linear term of this expansion is calculated. It is also shown that reduced distribution 
functions of order less than s don't possess such property and, to all appearance, 
contain not all information about the system under consideration. 
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1. Introduction 

In classical statistical mechanics an equilibrium system of molecules in a volume V 
is described by canonical distribution function Ff^{q,p), where is a set of phase 

variables: coordinates qi and momenta pi of molecules. If interaction of molecules 
is additive, reduced distribution functions are introduced [H E]. They are used for 
evaluation of thermodynamic characteristics of such molecular system. It's usually 
accepted that reduced distribution functions contain information about a molecular 
system less than the canonical Gibbs distribution function. It's also supposed that the 
lower an order of a reduced distribution function is, the less information it contains. 
But there does not exist a proof of this statement in scientific literature. 

On the other hand, it is known that for an equilibrium canonical ensemble of non- 
interacting particles a canonical distribution function F/v(g,p) is decomposed into a 
product of reduced one-particle distribution functions Fi{q,p) [IJ. This means that 
all information about such system is contained in the reduced one-particle distribution 
function. In [3l H] it was proved that for a closed molecular system with pair interaction 
there is a one-to-one correspondence between a canonical distribution and a reduced 
two-particle distribution function. 

In this paper it's proved that for a system having interaction potentials up to order 
s there exists one-to-one correspondence between full canonical distribution function 
and a reduced s-particle distribution function. This means that the s-particle function 
contains the whole information about system under consideration. 

We consider an equilibrium system of particles contained in the volume V under 
the temperature T. Potential energy of system is supposed to have the form 

s 

UN{qi,...,qN) = Y1 Yl M<lji,---,Qji), (1) 

1=1 l<ji<-<ji<N 

where s is an arbitrary fixed integer less than and . . . , qi) is a direct interaction 

energy of / particles. Probability distribution function of equilibrium system is the 
canonical Gibbs distribution which is decomposed into a product of a momentum 
distribution function and a configurational one [H [2]. The former is expressed as a 
product of one-particle Maxwell distributions, the latter has the form 

Djsfiqi, . . . , g^) = Q-^ exp{-/3f/7v(gi, . . . , gTv)}, (2) 

where f3 = 1/kT, k is the Boltzmann constant and Qa? is the configuration integral 

Qn = jexp{-pUN{qi, . . .,qN)}dqi ■ --dqN- (3) 

Here and below integrating with respect to every configurational variable is carried out 
over the volume V. For a system having interaction of form ([T]) reduced distribution 
functions are introduced by expressions [Ij 

A^! f 

Fi{qi, ...,qi)= _ J DN{qi, . • . , qN)dqi+i ■ ■ ■ dqN, / = 1, 2, .... (4) 



On completeness of description of an equilibrium canonical ensemble 



3 



These functions are used instead of full canonical distribution ([2]) to calculate various 
characteristics of the molecular system. Let us investigate properties of the reduced 
s-particle distribution function. 

Potential energy ([1]) can be written as 



l<ji<-<js<N 



where 



. . . , g,) = ^(C^_;.) ^ Yl ^^(^ii' ■ • ■ ' ^^r) (6) 

r=l l<ii<--<jr<s 

and = m!/{n!(m — n)\} are binomial coefficients. 
Introduce a function h{qi, . . . ,qs)} by the relation 

exp{-/30(gi, . . . , gj} = cr{l + . . . , g^)}, (7) 

where 

^ y exp{-/30(gi,...,g,)}rfgi---cig,. (8) 
The canonical Gibbs distribution (j2]) takes the form 

DN{qi, ...,qN) = Q~N^ n + ^'^^Ji' • • • ' 

i<ii<-<is<iv 

with 

Qn = Yl + ^'^^Ji' • • • ' ■ ■ ■ dqN- (10) 

l<ii<-<js<Af 

From and (fTOl) it follows that statistical properties of system under consideration 
are completely determined by the specifying single function of s configurational variables 
h{qi, . . . ,qs) and two external parameters N and V. The assumption naturally arises 
that another single function of s configurational variables -F5(gi, . . . , g^) can completely 
determine all statistical properties of this system. It turns out that this is indeed the 
case and reduced distribution function of order s plays a special role among all reduced 
functions (jlj). It can be proved that canonical Gibbs distribution ([2]) is expressed in 
terms of Fs(gi, . . . , g^). This means that there is one-to-one correspondence between 
D]\r and Fg. Therefore the system under consideration can be completely described by 
both canonical distribution and reduced s-particle distribution function. All functions 
(jlj) are expressed in terms of Fs(gi, . . . , g^) too. We prove these statements below. That 
proof is analogous to one for system with two-body interaction (s = 2) presented in the 
papers [HI Hj . 

In section |2] we pose a mathematical problem for our molecular system and 
formulate conditions for existence and uniqueness of its solution. In section [3] feasibility 
of these conditions for considered physical system are proved. In section Ulan expression 
for function /i(gi, . . . , g^) in terms of F5(gi, . . . , g^) is calculated. In section [5] an 
expression for the canonical distribution in terms of reduced s-particle distribution 
function is produced. In section [6] it is shown that reduced distribution functions of 
orders less then s don't possess such property. 
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2. Mathematical formulation of problem 

Let us introduce a function /(gi, . . . , g^) by the relation 

iV' 

^^(^1' • • • ' = (ATT^yr^ti + /(^i' • • • ' ^^)]- (11) 

Both the function /(gi, . . . , g^) and the function h{qi, . . . , g^) satisfy the conditions 

j /(gi, ...,qs)dqi--- dq^ = 0, ^ /i(gi, . . . , gj rfgi ■ ■ ■ dg^ = 0. (12) 
From expressions ([9]), and (fTT]) it follows that 

l + /(gi,---,gs) = 7— / II [l + /i(gji,...,gjj]c?g,+i---rfg^.(13) 

i<n<-<]s<N 

This relation defines the transformation {h ^ f} and can be considered as a nonlinear 
equation relative to /i(gi, . . . , gs). If there exists a solution /i(gi, . . . , g^; [/]) of this 
equation then the function D^iqi, ■ • • , gA?) becomes an operator function of /. It means 
that both the canonical Gibbs distribution and all reduced distribution functions 
Fi are expressed in terms of the single reduced distribution function Fs. Thus we have 
to prove that equation (fT3|) has a unique solution and therefore the transformation 
{h /} has inverse one {/ h}. 

Multiplying equation (fT3l) by QnV"^ and using (|TOl) we rewrite it in the form 

[^ + fiQi,---,Qs)]y]^ I n {^ + h{Pji,---,Pjs)]dpi---dpN 

l<ji<-<js<N 

-yW^ / n [1 + ^(^-Ji' • • • ' ^js)Ws+i ■■■dqN = 0. (14) 

l<ji<-<i,<iV 

The left-hand side of (1141) is a polynomial operator of degree N" = Cf^ relative to h and 
of degree one relative to /. Denote this operator by T{h,f). Equation f|T^ can be 
written in a symbolic form 

m/) = 0. (15) 

To solve this equation it is necessary to specify an additional condition 

jr(/,(o)j(o)) ^0, (16) 

where /'•°^(gi, . . . , g^) and h^^\qi, . . . ,qs) are assigned functions. 

We can easily determine these functions for our physical system. If the external 
field and all interactions between particles are absent, i.e., all potentials ui are constant, 
then the function /i(gi, . . . , g^) vanishes. Under this condition = Q^^^ = V^, 
D^lqi, . . . , q^) = D^^^ = , and /(gi, . . . , g^) = 0. Therefore we can take h^^'' = 
and /(°) = in (fT6|) 

Equation f|T5|) and additional condition (|T6l) form a problem on implicit function. 
In functional analysis there is a number of theorems on implicit function for operators of 
various smoothness classes. We use the theorem for analytic operator in Banach space 
in the form given in the book [5] 
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Theorem. (On implicit function). Let T{h,f) be an analytic operator in 
Dr{h^^\Ei) X Dp{f^^\E) with values m E2. Let an operator B = -dJ^{h'-°\ f^°^)/dh 
have a bounded inverse one. Then there are positive numbers pi and ri such that a 
unique solution h = xif) of the equation J^{h, /) = with the additional condition 
!F{h'^^\ /*^°^) = exists in a solid sphere Dr^{h^'^\ Ei). This solution is defined in a solid 
sphere Dp-^(f^'^\E), is analytic there, and satisfies the condition h^^^ = xif^^^)- 

Here Dr{xo,S) denotes a solid sphere of radius r in a neighborhood of the element xq 
in a normalized space S, the symbol x denotes the Cartesian product of sets, dJ-'/dh 
is Frechet derivative [3, [6] of the operator J-', h^^^ and /^°^ are assigned elements of the 
respective spaces Ei and E. If the functions h{qi, ■ ■ ■ ,qs) and /(gi, ■ ■ ■ ,qs) satisfy this 
theorem conditions, the former is a single valued operator function of the later. 

To prove an existence and uniqueness of a solution of problem (fT5l) . (JT6l) we have 
to show that all conditions of the above theorem are satisfied. 

3. Proof of feasibility of the theorem conditions 

First we define spaces Ei, E, and E2 mentioned in the theorem for functions describing 
the physical system under consideration. 

3.1. Functional spaces of problem 

Potentials Uk{qi, ■ ■ ■ , qu) are real symmetric functions. Suppose they are bounded below 
for almost all {gi, . . . , qu} G V . All physically significant potentials possess this property. 
Under this condition integral ([8]) exists and ... ,5^) is a real symmetric function 
bounded for almost all {gi, . . . , g^} G V . 

A set of functions bounded nearly everywhere forms a complete linear normalized 
space (Banach space) with respect to the norm |6l [7j 

\\h\\= vraisup . . . , g^)] , (17) 

{qi,...,qs}&V 

where "vrai sup" denotes an essential upper bound of the function on the indicated set. 
It is called the space of essentially bounded functions and is denoted by Loo(V''*-'). Here 

=^ X ■ ■ ■ X is a repeated s times Cartesian product of V by itself. In addition 
/i(gi, . . . ,qs) satisfies condition (fT2l) . The set of such functions is a subspace of Loo(V*-*^). 
It is easy to show that this subspace is a complete space relative to norm ( fTTl) . Therefore 
we can take the Banach space of symmetric essentially bounded functions satisfying 
condition f|T2l) as Ei. 

Expression ( |T3i) for /(gi,...,gs) includes multiple integrals of different power 
combinations of h{qi, . . . , g^). Any power of essentially bounded function are integrable 
with respect to arbitrary set of variables {g^j, . . . ,qj^} over V [8]. Therefore all integrals 
in f|T^ are essentially bounded functions too. Arguing as above, we can show that the 
space E of functions /(gi, • • • ,gs) coincides with Ei. Continuing in the same way we 
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can show that E2 is the same space. Thus we define the spaces of the above theorem as 
E = Ei = E2 = Loo(V(")) with property (fT2|). 

From (I7j) and ( ITTi) for h{qi, ■ ■ ■ ,qs) and /(gi, ■ ■ ■ ,qs) it follows that / > — 1 and 
h > —1. Therefore we can take a manifold {/>— 1, /i>— l}asa definition domain 
of the operator J-'{h, /). Since the left-hand side of ( IT4l) is a polynomial, the operator 
!F{h, f) is analytical in this domain. As stated above the additional condition f|T6|) is 
valid for h^^^ = f^^^ = 0. Thus any solid spheres of -Ei, E with centers at h^^^ = 0, 
Z'^") = and radii r < 1, p < 1 respectively can be used as domains Dr{h^^\ Ei) and 
Dp{f^^\ E) indicated in the theorem. 

Finally it is necessary to prove that the operator 
dof dJ^{h, f) 



B 



dh 

has a bounded inverse one. 



(18) 



3.2. Properties of the operator B 

To find the inverse operator B^^ it is necessary to solve the equation Bh = y, where 
h G El, y G E2. The expression for Bh is a linear relative to h part in the left-hand side 
of relation (IT^ as / = 0. Expanding products in (IT^ and keeping linear summands we 
obtain 

{Bh){qi,...,q,)) = h{qj^,...,qjJdqs+i---dqN. (19) 

i<ji<-<js<N 

Variables {g^^, ■ ■ ■ ,qjs} are divided into two groups {gj^, . . . , gj,} and {gj,+i, • • • , gj^} fo'^ 
/ < s, where (ji, ■ ■ ■ ,ji) C (1, . . . , s) and {ji+i, . . . ,js) C (s + 1, . . . , N). The sum in the 
right-hand side of (fT9l) is divided into two parts for each < / < s. The expression (fT9ll 
can be rewritten in the form 

{Bh){qi,...,q,)=j2 E E 

i=0 l<ji< - <ji<s s+l<ii+i<-<is<Af 

■-^7^ y h{qj„ . . . ,qj„qj,^^, . . . ,qjjdqs+i - ■ ■ dqN- (20) 

We suppose that s < A^/2. Integration in ( l20i) is carried out with respect to 
variables from second group {g^j^^, . . . , g^^} C {g^+i, • • • , gA^} and remaining ones 
{qs+i, . . . ,qN}\{qji+i, ■ ■ ■ ,qjs}- Here a symbol \ denotes difference of sets. Since 
the function h{qi, . . . ,qs) is symmetric, we see that the summation with respect to 
{ji+i, . . . ,js) gives C^l'^ identical terms. In the result we get 

s 

{Bh){qu...,qs)=J2C'NJs E 

1=1 l<ji<-<ji<s 

■-^TJry j h{qj^, . . . ,qj^,pi, . . . ,ps-i)dpi ■ ■ ■ dps-i. (21) 
Here we omitted summands with / = because of property (fT2l) for the function h. 
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It's easy to estimate a norm of the operator B. Using definition (fTTI) we obtain 

s 

\\Bh\\<Y,C^^-J, \\h\\ 

1=1 l<ji<-<ji<s 
s 

= J2^N-sCi\\h\\ = (C^ - C;,_Mh\\. (22) 
(=1 

From here we get an estimation 

||5||<(C^-C^_J. (23) 

Thus the operator B is bounded. 

To simphfy transformations we introduce notation 

h (gi, . . . , gfc) = / h{qi, qs)dqk+i ■ ■ ■ dq^. (24) 

Using ( l2Ti) and (12^ we can write the equation Bh = y in the next form 

s 

Y.^N^ E h^\qn,...,qj,) = y{qu...,qs). (25) 

1=1 l<ji<-<ji<s 

Integrating this equation wit respect to g^, then qs_i and so on we obtain a system of 
equations for h^''^ 

k 

Y.^N-k E h^\qn,...,qj,)=y^'\qi,...,qk) (26) 

1=1 l<ji<-<ji<k 

for all/c = 1, . . . , s. Correctness of this expression for arbitrary k is easily tested by 
induction. Note that = and h^'^^ = h. 



A solution of system (|26|) is obtained step by step starting from the first equation 
with k = 1. It has the form 

h"'\qi, . . . , g.) = E(-l)'=-^^y^^ E y^'^i^n, • • • ' ^> )(27) 

for all = 1, . . . , s. This solution is easily checked by forward substitution into system 



Putting here k = s we obtain the expression for the inverse operator B ^ 
h{qi, ...,q,) = {B'^y){qi, ...,qs) 



r=l l<ji<---<jr<s 

From here it's easy to estimate the norm of the inverse operator B~^. Evaluate the 
norm of right-hand-side of ( l28i) 



N-r ^ " " ^N-r 

T=l l<ii<--<ir<s r=l 
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yN- 



We obtain from here 

^ A/" — s 

\\B'^\\<yT^C:<2'-l. (30) 

r=l 

Therefore the operator exists and it is bounded. 

So all conditions of the above theorem are valid for our physical system. Hence 
there exists a unique solution h of problem ([T5|) . f|T6l) as a function of /. This 
solution defines an inverse transformation from the function / to the function h: 
h{qi,...,qs) = h{qi,...,qs] [/]). 

4. Derivation of the inverse transformation h{f) 

To obtain the transformation f h we have to solve equation ( fT4|) relative to 
h{qi, . . . ,qs). At first define an auxiliary operator function g{h) by means of a relation 

^ n [^ + h{(lji,---,(ljs)]dQs+i---dqN = l + g{h). (31) 

l<ji<-<js<N 

This operator function is a polynomial of degree A/" = relative to h and depends on 
s configurational variables {gi, . . . , g^}. It can be written in the form 

g{h)=J29lih), (32) 

1=1 

where gi{h) is a uniform operator of order / relative to h. Let us derive the expression 
for gi{h) from fl3Tl) . 

For the sake of abbreviation of subsequent calculations we introduce next notations. 
We will denote by number K G {l,...,J\f) every ordered collection (ji,...,js) C 
(1,...,A^). Such one-to-one correspondence can be always made. A collection 
{qj-^, . . . ,qjj is an element of manifold V'-'^-*. We will denote this element by Xk- By 
definition put Xi = (gi, . . . , g^). 

Expanding the product in ( l3Ti) we obtain for every I = 1, . . . ,J\f 

9i{h) = — ^ f dqs+i ■■■dqN ^ KXk^) ■ ■ ■ h{XK,). (33) 

l<Ki<-<Ki<M 

Introduced operators g{h) and gi{h) are symmetrical functions of s configurational 
variables: g{h) = g{qi, . . . ,qs][h]) and gi{h) = gi{qi, . . . ,qs][h]). In contrast to h 
and / both g{qi, . . . ,qs;[h]) and (//(gi, . . . , g^; [h]) don't satisfy condition f|T2|) except 
for (7i(gi, . . . , g^; [/i]). First term of series fl52]) is (7i(/i) = -B/i and satisfies to condition 
( IT2|) . Configuration integral (ITOl) takes the form 

Q^ = r^[l+^W]=l^^[l + f;^(°)]. (34) 

fc=2 

Here we used notation ( l24l) . 

Substituting definitions (!3T]) and into (IT^ we write it in the form 

(l + /)[l + ^W(/i)]-[l + ^?(/i)]=0, (35) 
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where / and g{h) are functions of Xi = (gi, . . . , g^)- But value g^^^ doesn't depend on 
configurational variables. It is a functional relative h. Substituting here expansions (l32l) 
and ( IMl) we reduce this equation to the form 

M M 

/[I + Y.9f\h)] + Y.[9f\h) - -Bh = 0. (36) 

1=2 1=2 

Here we have taken into account that gi{h) = Bh and 'g^^\h) = 0. 
For subsequent calculation we need multihnear operators 

Gi{yi,...,yi) = I dqs+i---dqN ■ ■ ■ y,(X;^J. (37) 

i<Ki<-<Ki<j\r 

These operators are linear with respect to any functional argument y,. We can consider 
the operator functions gi{h) as generated by these multilinear operators Gi 

giih)=Gi{h,...,h). (38) 

Operators Gi{yi, . . . ,yi) are functions of configurational variables (gi, . . . , g^) = ^i- 
In general these functions aren't symmetrical relative to (gi,...,gs). But this isn't 
important since under substituting of these operator functions into equation (l36l) 
symmetric property will be hold automatically. In the result we can rewrite equation 
§6^ as 

h = B''f[l + Y,Gf\K...M 

1=2 

N 

+B~' ^[Gr(/^, . . . , /i) - Giih, . . . , h)]. (39) 

1=2 

We will search a solution of this equation in the form of power series 

oo 

h = Y,h{f), (40) 
fc=i 

where hk{f) are uniform operators of order k relative to /. At the same time they are 
functions of configurational variables Xi. Substituting fHOj) into fl39l) and taken into 
account linearity of Gi{yi, . . . , yi) with respect to any argument yi we obtain 

oo A/" oo oo 

huif) = B-'f + E E ■ ■ ■Y.^i'^hn, • • • , h,,)f 

k=l 1=2 ji=l ji=l 

M oo oo 

E E • • • Y}^f\hn^ ■■■,h,,)- Giih,,, . . . , h,,)]. (41) 

1=2 ji=l 31=1 

Transform sums over ji , . . . , ji as follows 

oo oo oo 
Ji=l ii=l k=l jiH Vii=k 
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Then relation (1411) takes the form 

oo N oo 

Y.hk{f) = B-'f + B-'Y.Y. E Gf\hn^....h,;)f 

k=l 1=2 k=l jiH \-ji=k 

M oo 

E [Gf\h,,,...,h,,)-Gi{h,,,...,h,,)]. (43) 

1=2 k=l ji + -+ji=k 

Double sum YJ1L2 ^"kLi transformed as follows 

Af 00 Af-lAf~l Af 00 AT-l k oo Af 

1=2 k=l 1=2 k=l 1=2 k=jV k=2 1=2 k=jV 1=2 

or as follows 

Af oo Af Af Af oo Af k oo Af 

EE = EE+E E =EE+ E E- 

1=2 k=l 1=2 k=l 1=2 k=Al+l k=2 1=2 k=Af+l 1=2 

Substituting (144 gp and (144 6p into (|43l) we get the relation 

oo Af fe— 1 

J2 hkif) = B-'f +B-'Y.Y. E Gf\h,,, h,^)f 

k=l k=3 1=2 jiH hji=A:-l 

oo Af 

E E E G^\h,,,...,h,)f 

k=Af+l 1=2 ji+--+ji=k-l 
Af k 

+5-' EE E [Gf{h,,,...,h,,)-Gi{h,,,...,h,,)] 

k=2 1=2 ji + -+ji=k 
oo Af 

+B-' E E E [Gf{h,,,...,h,^)-G,{h,„...,h,^)]. (45) 

k=Af+l 1=2 ji+--+ji=k 

Here in the first two sums of the right-hand side we change summation variable k to 
k + 1. 

In this relation all sums with respect to k contain expressions of order k relative 
to /. Putting terms of the same order being equal in accordance with the theorem 
on uniqueness of analytical operators [9] we obtain the next recurrent system for the 
functions hk{f) 

hi{f)=B-'f, (46a) 
h2if) = B-'[Gf\h,h) - G2ih,h)], (466) 

fe-l k 

h{f) = B-'Y. J2 Gf{h,„...,h,,)f+B-'Y1 E 

1=2 ji+-+ji=k-l 1=2 ji+-+ji=k 

■[Gl'\h^,, . . . , /i, J - Giih^,, . . . , h,^)], 3<k<Af, (46c) 

Af Af 

hk{f) = B-'Yl E G^'\hn,---,h,,)f + B-'Yl E 

1=2 ji+-+ji=k-l 1=2 ji+-+ji=k 

■[Gf{h,,, . . . , /I, J - Gi{h,,, . . . , k>Af+l. (46d) 
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All terms of series f l4UI) are calculated from this system. So the solution of the equation 
f|T5|) is founded. It satisfies additional condition f|T6|) . Convergence of series fHOl) with 
hk{f) being the solutions of system f l46al) -( l46dj) is proved by Cauchy-Goursat method 
presented in the book [5]. 

As soon as hk are expressed in terms of / we can get the canonical distribution ([9]) 
in terms of Fg since / and Fg are uniquely connected by relation (fTTj) . 



5. Calculation procedure for the canonical Gibbs distribution 

Since canonical distribution ([9]) is a ratio of two polynomials with respect to h, we see 
that Dn is an analytical operator function of h. We have just proved that h is an 
analytical operator function of /. Therefore Dj^ is an analytical operator function of / 
and it can be expanded into an absolutely convergent series relative to / 



l + J2Mf)]^ (47) 



k=l 



where fkif) is a uniform operator of order k transforming function f{qi,---,qs) to 
function fkili, ■ ■ ■ , Qn] [/])• Taking into account the definitions of reduced distribution 
functions (jlj) and function /(gi, ■ ■ ■ ,qs) ( pTl) we can get relations for (pk 

dqs+i---dqN<^i{qi,...,qN;[f]) = f{qi,---,qs), (4^ 



yN- 

dqs+i- ■ ■ dqN^Pkiqi, ■ ■ ■ ,qN][f]) = 0, k = 2,3,.... (49) 



Below we construct a procedure for calculation of functions ipi{qi, . . . ,qN] [/]) in terms 
of/. 

We introduce a nonlinear operator function X{h) by the relation 

n [l + %,„...,g,J] = l + A(/i). (50) 

l<ji<-<js<N 

This operator function is a polynomial of degree Af relative to h. It can be written in 
the form 

M 

A(/i) = ^A,(/i), (51) 

k=l 

where Xk{h) are defined by relations 

\kiqu...,qN;[h])= Yl h{XK,)---h{XK,). (52) 

l<Ki<-<Kk<N' 

Introduce also multilinear operators 

Afc(yi, . . . , 2/fc) =^ 5Z Vii^Ki) ■ ■ ■ VkiXxJ- (53) 

l<Ki<--<K^,<J\f 

It's evident that 

Xkih)=Akih,...,h). (54) 
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The operators introduced here are connected with the operators g{h), Qkih), and 
Gk{hi, . . . , hk) by the relations 



yN 



yN 
1 



j dqs+i ■ ■ ■ dqNXiqi, • • • , g^v; [h]) = g{qi, . . . , g,; [/i]), (55) 

rfg^+i ■ ■ ■ dq^Xkiqi, ...,qN][h])= gk{qi, ■ ■ ■ , g^; [h]), (56) 



dq,+i ■ ■ ■ dq^Akiqi, . . . , giv; [h, ...,hk]) 



yN- 



yN- 

= Gkiqi,...,qs;[hi,...,hk\). (57) 
In particular for A; = 1 

dqs+i ■ ■ ■rfg7vAi(gi, ...,qN; [h]) 

= Gi{qi,...,qs;[h]) = ^i(gi, . . . , g,; [/i]) = (B/i)(gi, . . . , g,). (58) 
Taking into account the expression flMl) for Qat we can write 

Comparing it with (H7|) we get the relation 

where /i is the operator function of / calculated in previous section. Using here the 
expressions for X{h), Ij^^^h) and h[f) we can transform the right-hand side of fl60l) 
to series with respect to / and thus obtain expressions for fkif)- But less awkward 
transformations are obtained if we construct a recurrent system for fkif)- 
Multiplying ([60]) by l + g^^\h) and using §Bl and (JSl) we obtain 

AT oo Af N 

fc=2 1=1 k=l k=2 

Substitution of the expansion (HDl) here gives 

J\f oo oo oo 

k=2ji=l jk=l 1=1 

AT oo oo AT oo oo 

= E E • • • E Mhn,-, - E E ■ ■ ■ E Gk\hn, /^.J- (62) 

k=l ji=l jk=l k=2 ii=l jk=l 

Further calculation is carried out in the same way as in the previous section. We 
won't make it and write a recurrent system for fkif) straight away 

<^i(/) = Ai(/ii), (63a) 

^2{f) = Ai(/i2) + A2{h, h) - (/ii, h), (636) 

k 

Mf) = Mhk) + Yl E {Hhn,-,hj,)-Gf{h,,,...,h,,)} 

1=2 ji + -+j,=k 
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k 

-E E G\'],{h,„...,h,,_,)y,,,{f), 3<k<Ar, (63c) 

'=3 jiH \-ji=k 

N 

Mf) = Ai(M +E E {Mhn,-.h)-Gfihn,-,h,,)} 

1=2 ji+-+ji=k 
Af+1 

-E E k>^^+l. (63c?) 

'=3 iiH hi;=A: 

In this relations we have to use the expressions for hr{f) derived from the recurrent 
system fH6al)- fl46dl) . 

For example an expression for v^i(/) is 

s 

^i(gi, . . . , g^; [/]) = E(-1)^''C'— „i E /^'^(^^- • • • ' (^4) 

r=l l<jl<--<ir<Af 

It is easy to show that expressions fl63fflp - fl63Gjj) and flMl) satisfy conditions fllS]) . fH^ . 
For s=2 expression coincides with (y9i(gi, . . . , g^r; [/]) derived in the papers [3l H]. 

6. Inadequacy of reduced distribution functions of order less then s 

Let us test the theorem conditions for reduced /-particle distribution function when 
/ < s. Introduce a function //(gi, . . . , qi) by relation 

Fi{qi, . . . , qi) = ^^_^^,^J l + /;(gi,--.,gO]- (65) 

All constructions and reasonings of sections[2]and[3]remain valid. We obtain the operator 
equation J^i{h, fi) = and appropriate additional condition. There are proper Banach 
spaces and bounded operator Bi which is Frechet derivative of the nonlinear operator 
J-i{h, fi) relative to h. It is easy to show that the operator Bi is resulted by integrating 
of the operator B with respect to g^+i, . ■ ■ ,qs- Therefore a uniform equation Bih = 
has the form (see derivation of equation fl2^ ) 
I 

Y.Cn~i E h^''\qn^---,q,r) = 0. (66) 

r=l l<ji<--<jr<l 

From here we see that the operator Bi has a nontrivial space of zeroes. This space 
consists of all functions h{qi, . . . , g^) G -Ei satisfying a condition 

J h{qi, qs)dqi+i ■ ■ ■ dqs = (67) 

for any fixed /. This means that the operator Bf^ isn't exist as / < s and inverse 
transformation {fi h} isn't exist either. Therefore we can't express the canonical 
distribution Z^a? in terms of as Z < s. As soon as / = s the operator Bi = B and the 
space of zeroes of operator Bi becomes trivial (see the condition fl^TI) ). In this case all 
conditions of the theorem are valid and we obtain all above results. So the reduced s- 
particle distribution function plays a specific role for the system with s-body interaction. 
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Fg is a reduced distribution function of minimal order containing all information about 
this system. 

7. Conclusion 

Using the theorem on implicit functions in this article it is shown that a reduced 
distribution function of order s plays a specific role for a canonical ensemble of N 
particles with s-body interaction. The canonical Gibbs distribution Dj^f{qi . . . , q^) 
can be expressed uniquely in terms of this function Fs{qi . . . , qs) ■ From here we 
easily conclude that there is a one-to-one correspondence between these two functions. 
This means that the reduced distribution function Fg contains information about 
system under consideration as much as the whole canonical distribution Djy. Reduced 
distribution functions of all orders can be expressed in terms of this single function F^. 

All reduced distribution functions of order I less then s don't satisfy the theorem 
conditions. So it is impossible to express the canonical distribution in terms of these 
functions of order I < s. To all appearance they contain not all information about the 
system under consideration. 

Considered theorem provides sufficient conditions for existence and uniqueness of 
inverse transformation {/ h}. Results obtained here are valid in some neighbourhood 
of /i(o) = 0, /(°) = 0. The question about size of this neighbourhood demands special 
investigation. 
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